Breakdown of half-metallic ferromagnetism in zinc-blende II-V compounds 
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We investigated the electronic and magnetic properties of a series of zinc -blend II-V compounds by carrying 
out density-functional-theory calculations including spin-orbit couplings. Contrary to the case of CaN and CaP, 
the half-metallic characteristics of the II-V compounds such as CaSb and CaBi were found to be destroyed. 
Our analysis of the valence band structures of CaAs, CaSb, and CaBi revealed a critical role of the spin-orbit 
coupling interactions on the exchange-split band structure, thereby leading to breakdown of the half-metallic 
ferromagnetism for the systems with heavier group V elements in the zinc-blend II-V compounds. 
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Half-metallic ferromagnets (HMFs), being metallic in only 
one of the two spin channels, have been considered as an in- 
dispensable ingredient in the development of spintronic de- 
vices and applications. The realization of the half-metallicity 
has been investigated extensively since its first prediction 
on Heusler alloys by de Groot et al.—. However, due to 
the complexity in their electronic characteristics, manifesting 
both metallic and insulating properties in a single system at 
the microscopic level, HMF materials often require complex 
structures such as ternary spinel, Heusler, and double per- 
ovskite structures, also including transition-metal elements as 
a source of local magnetic moments. 

To explore possible spintronics applications of HMFs to 
semiconductor devices, there have been a great deal of stud- 
ies exploiting the half metallicity in the binary compounds of 
zinc-blende (ZB) structur e) 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 , which is sim- 
ple and compatible with existing III-V and II- VI semiconduc- 
tors. So far there have been a few reports on the fabrication of 
nano-scale MnAs dots on GaAs, CrSb ultra-thin films on the 
GaSb substrate, and ultra-thin CrAs layers in the CrAs/GaAs 
multilayer s - 2,3 ' 4 ' 5 . In addition, the half-metallic ZB com- 
pounds with transition-metal elements were also reported by 
first-principles density-functional-theory (DFT) calculations, 8 
whereas a high Curie temperature above 400 K has been ob- 
served in experiments for CrAs and CrSb systems 2,3 . 

Contrary to the ZB compounds containing transition-metal 
elements, where the localized d electrons are responsible for 
the ferromagnetic component in HMF, several DFT calcula- 
tions have predicted another kind of ZB II-V compounds as 
a candidate for the HMFs without containing any transition- 
metal elemen t 10 ' 12,13 . The ZB compounds of alkaline earth 
elements Ca, Sr, and Ba combined with all the elements of 
the group V were shown to be HMFs with a magnetic mo- 
ment of one Bohr magneton (/j,b) per formula unit-cell (f.u.). 
In these compounds, the presence of a flat p band crossing the 
Fermi level (Ep) in its paramagnetic phase is a key to the half- 
metallic electronic structure. The narrow p band contributes 
to the exchange energy splitting close to 0.5 eV, consequently 
leading to an insulating gap in the majority spin channel and 
a metallic state in the minority spin channel. 

However, since the p states of the heavy elements like 
Sb and Bi are affected by the relativistic spin-orbit coupling 
(SOC), the SOC is expected to play a role in the determina- 
tion of the II-V valence band structure consisting of the an- 



ion p components of As, Sb, and Bi atoms. Therefore, their 
ground state electronic and magnetic properties needs to be 
examined carefully in connection with the half-metallicity of 
the II-V compounds with heavy group V elements. 

In this paper, we report the results of our non-collinear DFT 
calculations including the spin-orbit coupling terms for the 
electronic and magnetic properties of the ZB II-V compounds 
of AX (A=Ca, Sr, and Ba; X=N, P, As, Sb, and Bi). Contrary 
to the previous reports on the ZB II-V compounds, we found 
that the ZB II-V compounds containing Sb or Bi are no more 
half-metallic due to the strong SOC overriding the exchange 
instability triggered by the flat band feature of the anion p 
states. Our calculations showed a broad spectrum in the phys- 
ical properties of the ZB II-V compounds, ranging from the 
half-metallic ferromagnet to an ordinary ferromagnetic metal 
to a paramagnetic metal. Further analysis of the valence band 
structures explains the competition between spin-orbit cou- 
plings and exchange interactions as a key ingredient in the 
determination of ground states of the ZB II-V compounds. 

To investigate the role of spin-orbit coupling interactions, 
we carried out non-collinear DFT calculations including the 
SOC to determine the electronic structure and magnetic prop- 
erties of the ZB II-V compounds of AX (A=Ca, Sr, and Ba; 
X=N, P, As, Sb, and Bi). For the self-consistent electronic 
structure calculations, we have optimized the equilibrium lat- 
tice constant of each AX compound under the symmetry con- 
straint of the zinc-blende (ZB) structure. The structural stabil- 
ity has been extensively examined in the previous worki^. It 
was discussed that the ZB-type structure of CaAs can exist as 
a metastable state under certain conditions even though nei- 
ther ZB nor wurtzite (WZ) structures are the absolute ground 
state of CaAs. Without further discussions on the stability of 
the ZB structure, from now on, we like to focus on the elec- 
tronic and magnetic properties of the II-V compounds at their 
equilibrium positions within the ZB structure. 

The electronic structure calculations were performed by 
using the OpenMX DFT code 14 , a linear-combination-of- 
pseudoatomic-orbitals method 15 . The generalized gradient 
approximation (GGA)i£ was used for exchange-correlation 
potential. The effects of spin-orbit couplings via a relativis- 
tic j-dependent pseudopotential scheme were included in the 
non-collinear DFT calculation s 17 ' 18 - 19 . Double-valence-plus- 
single-polarization orbitals were taken as a basis set, which 
were generated with cutoff radii of 6.0, 6.0, and 7.0 a.u. for 
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(a) "SOC off" CaN 



(b) "SOC on" CaN 



FIG. 1: Total density-of-states of CaX (X = N and P) without and 
with spin-orbital coupling, as labeled by "SOC off" and "SOC on", 
respectively. 



FIG. 2: Band structures of CaN (a) without ("SOC off") and (b) with 
("SOC on") spin-orbital coupling. 



Ca, Sr, Ba atoms, and 5.0, 5.5, 6.0, 6.5, and 7.0 a.u. for 
N, P, As, Sb, and Bi atoms, respectively. Troullier-Martins- 
type pseudopotentials, with a partial core correction for all 
atoms, were used to replace the deep core potentials by norm- 
conserving soft potentials in a factorized separable form with 
multiple projectors. The real space grid techniques were used 
with the energy cutoff of 300 Ry in numerical integrations 
and the solution of the Poisson equation using fast Fourier 
transformations. For the formula cell containing two atoms, 
(16x16x16) k-grid was sampled over the full Brillouin zone. 
To compare the stability of different magnetic configurations, 
we calculated total energies of paramagnetic (PM), ferromag- 
netic (FM), and anti-ferromagnetic (AFM) configurations in a 
supercell containing 16 atoms with (8x8x8) k-grid. 

Figure Q] shows the total density-of-states (DOSs) of ZB 
CaN and CaP with and without spin-orbit coupling (SOC). 
Hereafter we denote the results with and without SOC by us- 
ing the labels "SOC on" and "SOC off", respectively. The 
"SOC off" and "SOC on" DOS's for both CaN and CaP are 
almost identical. It is obvious that the SOC effect is negligi- 
ble for the light elements like N or P atoms. Consequently 
the HMF character with a magnetic moment of 1 /^/f.u. re- 
mains unchanged regardless of the presence of SOC. In fact, 
the spin-orbit coupling has virtually no effect on the ZB CaN 
band structure, as shown in Fig. [2] While the spin-up and 
spin-down bands are manifest in the collinear, i.e., "SOC off" 
band structure of Fig. a), the non-collinear band structure 
as obtained from the "SOC on" calculations requires a spin- 
resolved representation of the bands. In the spin-resolved 
band plots, as illustrated in Fig.|2jb) and other "SOC on" band 
plots in the following figures, we marked each band state by 
the symbols (A or y), the size of which corresponds to the 
weight of each spin-component. Figure |2jb) clearly demon- 
strates that the spin-resolved band structure of CaN with neg- 
ligible SOC has a clear separation of spin-up and spin-down 
bands, which is virtually identical to the collinear bands of 
Fig. Ha). 
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(a) "SOC off" CaAs 
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(b) "SOC on" CaAs 



FIG. 3: "SOC off" and "SOC on" band structures of CaAs. 



To demonstrate the effect of the SOC interactions on the 
valence band structures of the ZB II-V compounds, we have 
calculated the band structures of CaX (X=As, Sb, and Bi) in 
addition to CaN and CaP. The band plots of CaAs, CaSb, and 
CaBi are shown in Figs. [3] E and [3] respectively. One of the 
common features of the "SOC off" band structures of CaX 
is a strong exchange-splitting of the flat p bands consisting 
of the anion p states. The origin of magnetism for the HMF 
CaX compounds was attributed to the magnetic instability of 
the narrow p bands at the Fermi level' 2 . In fact, the flatness 
of the p bands are more pronounced in the "SOC off" band 
structures of CaX (X=As, Sb, and Bi) than those of CaN and 
CaP. 

The nature of the flat p bands in CaAs has been discussed in 
the previous work-abased on the tight-binding model. We can 
understand the valence band structure of AX by a simple fee 
lattice of anion X (X= N, P, As, Sb, and Bi) with an extraor- 
dinarily large lattice constant. In this simplified picture, for 
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FIG. 4: "SOC off" and "SOC on" band structures of CaSb 
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(a) "SOC off" CaBi (b) "SOC on" CaBi 



FIG. 5: "SOC off" and "SOC on" band structures of CaBi. 



instance, the valence band structure of AX can be mimicked 
by the partially filled p bands of a fee Br solid, which features 
the same p band with the same electron filling as that of AX 
except the flat bands observed in AX. Indeed, to acquire the 
flatness, the X p-p hybridization should be counter-balanced 
by the A d-X p hybridization. 

Figure|6]illustrates charge density profiles of the "SOC off" 
spin-up states of CaAs, which are labelled by I\, Li, and Xi 
(i=l,2,3) in Fig. [3ja). At the T point, the triply degenerate 
Ti (i=l,2,3) states consist of p x , p y , andp z orbitals and form 
a anti-bonding configuration among the As p states. Due to 
the presence of Ca, however, the p-p anti-bonding states are 
mixed with the Ca d-As, p bonding configurations. Without 
the Ca d-As p bonding contribution, the energy level of the Tj 
states should have lied at the higher position than those shown 
in Fig. [3ja). As the k vector moves along TX and TL, the 
1(1,2 and Xl,2 states change into a non-bonding configuration, 
while the L3 and A3 states develop into a strong As p-As p 
bonding configuration, respectively, as illustrated in Fig.|6jc) 



b)L 2 Mr, (h )x a(# ^ 



(C)L 3 






FIG 6: (Color online) Charge density profiles of (a) Li, (b) L2, (c) 
L 3 , (d) Ti, (e) T 2 , (f) F 3 , (g) Xi, (h) X 2 , and (i) X 3 states of the 
"SOC off" CaAs bands, which are labelled in Fig.Oa). Light (red) 
spheres represent for Ca atoms and blue (dark) spheres for As atoms. 



and (i). Since the nearly flat bands of the doubly degenerate 
1(1,2 — 1\ 2 — --Xi.,2 states contribute to the high DOS at Ep, 
when the SOC is turned off, the Stoner instability would drive 
the system into the ferromagnetic ground state. 

Now let us consider the effect of spin-orbit coupling inter- 
actions on the valence band structure. When SOC is turned on, 
contrary to the robust HMF band structures of CaN and CaP, 
the electronic structures of CaA (A=As, Sb, and Bi) change 
remarkably. Despite that the "SOC off" valence bands of all 
the CaA (A=As, Sb, and Bi) compounds are almost identical 
as shown in Fig.[3]-[5] the band dispersion and spin characters 
of the "SOC on" band structures are markedly different and 
evolve as the atomic number of the element X increases from 
As to Sb to Bi. Since the SOC strength becomes larger for 
the heavier atoms, the spin-resolved band structures exhibit 
more complicated features. Before discussing the details of 
the SOC effects on the spin-resolved band structures, let us 
examine the evolution of the "SOC on" DOS of CaA (X =As, 
Sb, and Bi). 

The "SOC off" DOSs of CaA (A=As, Sb, and Bi) in 
Fig- L7I a X ( c ) an d ( e ) nave basically the same characteristics 
as those of CaN and CaP in Fig. [TJa) and (c) except the pro- 
nounced flatness in the CaA DOSs. As the SOC strength in- 
creases from CaAs to CaSb, the sharp peaks corresponding 
to the flat p bands broaden and the spin-up and down compo- 
nents seem to smear into each other. (See Fig. 13d) of CaSb.) 
For the CaBi case, where the SOC energy scale is dominant, 
the DOS feature shows no spin-polarization but a clear sepa- 
ration of two peaks: one centered at —2.2 eV corresponding 
to the j=l/2 state and the other at —0.5 eV to ,7=3/2, as shown 
in Fig. Etf) of CaBi. 

In addition to the evolution of the DOS features, one can 
observe the change of spin-polarization (SP) at Ep as well as 
the magnetic moment (A/). The change of SP can be easily 
understood from the "SOC on" DOS plots of Fig. [7] The spin- 
polarization of CaAs is close to but not exactly 100%, which 
means that CaAs is not a HMF in a strict sense, whereas the 
spin-polarization of CaBi is zero. While the half-metallicity 
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TABLE II: Exchange and spin-orbit splitting energies in AX (^4=Ca, 
Sr, and Ba; X=As, Sb, and Bi). Energies are in eV unit. 
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FIG. 7: Total density-of-states of CaX (X=As, Sb, and Bi) without 
("SOC off") and with ("SOC on") spin-orbital coupling. 



TABLE I: Spin-polarization (SP) at Fermi level and magnetic mo- 
ment (M) ^i_g/f.u. as obtained from the "SOC on" calculations. 



N P As 

SP I M SP/M SP / M 



Sb Bi 

SP/M SP/M 



Ca 1.00/1.00 1.00/1.00 0.96/1.10 0.81/1.11 0/0 
Sr 1.00/1.00 1.00/1.00 0.97/1.10 0.84/1.17 0/0 
Ba 1.00 / 1.00 1.00 / 1.00 0.97 / 1.03 0.79 /LOO 0/0 



in CaAs and CaSb is destroyed moderately, the half-metallic 
character of CaBi disappears completely. 

Table U summarizes the spin polarization (SP) at Ep and 
the magnetic moment (M) of AX (A=Ca, Sr, and Ba; X=N, 
P, As, Sb, and Bi). Here the SP is defined by the ratio of 
the spin-up and spin-down DOS at Ep over the total DOS 
at Ep. The magnetic moments in Table J] include the contri- 
bution from both the spin and orbital moments as obtained 
from the non-collinear calculations. Due to the smearing of 
the spin-up and down components near Ep, both As and Sb 
compounds have non-negligible contribution of the majority 
spin DOS at Ep. The mixed spin-components in the valence 
bands are clearly shown in the spin-resolved band structures 
of Fig. |3jb) and Fig. |4jb). As the SOC strength increases, 
some of the non-collinear band states have a strong mixture 
of spin-up and down components. 

The "SOC on" band structures of CaX (X=N, P, As, Sb, 
and Bi) shown in Fig. [2}j5] can be classified into three cate- 
gories: (i) the CaN-type where the spin-exchange energy of 
AE ex w 0.5 eV is dominant, (ii) the CaBi-type where j=l/2 
and j=3/2 bands are well separated by the SOC energy scale 
of AE so w 2 eV, and (iii) the CaSb-type, an intermediate case 
where both AE ex and AE so compete with each other. In the 
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case of (i) the CaN-type, the SOC interactions are negligible, 
i.e., AE ex ^ AE so , so that the ferromagnetic band struc- 
ture is set by the exchange instability triggered by the narrow 
p bands. On the other hands, to understand the band struc- 
tures of the types (ii) and (iii), one has to consider both SOC 
and exchange interactions on the same basis. To estimate the 
magnitudes of AE ex and AE so , we calculated the exchange 
and spin-orbit splitting energies at the T-point for the "SOC 
off" ferromagnetic and "SOC on" paramagnetic states, re- 
spectively. The calculated values of AE ex and AE so for AX 
(A=Ca, Sr, and Ba; X=As, Sb, and Bi) are listed in Table HH 

For the case of (ii) the CaBi-type, the SOC energy of 
AE so = 1.88 eV is much larger than the effective exchange 
energy of AE ex = 0.41 eV. Since the SOC term dominates, 
the valence band eigenstates of CaBi at the T-point are well 
represented by j=3/2 and j=l/2 states, corresponding to the 
4-fold and 2-fold degeneracy, respectively. Even away from 
the T-point, the spin-up and down components are mixed up 
by maintaining its character of the spin-orbit coupled j-states. 
Consequently the average spin-polarization in CaBi becomes 
null so that the ground state is non-magnetic. 

The valence band structures of CaSb, SrSb, and BaSb are 
much more complicated as the energy scales of the SOC and 
exchange interactions become comparable to each other. For 
an example, as listed in Table HH the SOC and exchange ener- 
gies for CaSb are AE so = 0.52 eV and AE ex = 0.44 eV, re- 
spectively. The overall shape of the "SOC on" valence bands 
of CaSb in Fig.Hfb) resembles those of the "SOC off" valence 
bands in Fig.Ua), whereas the degeneracy of the bands along 
L — r — X is significantly destroyed. Due to the strong con- 
tribution of AE so , the spin-flip term (L + S- + L-S+) in the 
SOC term of H so = A so L • S acts as a perturbation to the 
exchange-split states. 

In conclusion, we have shown that the half-metallic ferro- 
magnetism in ZB II-V compounds can be destroyed for AX 
(A=Ca, Sr, and Ba; X=As, Sb, and Bi) due to the presence 
of strong spin-orbit coupling interactions. Our non-collinear 
DFT calculations revealed the evolving features of the spin- 
orbit coupling and exchange-split in the spin-resolved band 
structures as the atomic number increases from As to Sb to 
Bi. Considering the variation of the spin polarizations of ZB 
II-V compounds, it is suggested that one can control the SOC 
strength to design either HMF or paramagnetic metal. This 
SOC-dependent feature may be useful in the spintronics ap- 
plication of HMF materials. 
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